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Abstract 

Inflation is terminated by a non-equilibrium process which finally leads to a thermal 
state. We study the onset of this transition in a class of hybrid infiation models. 
The exponential growth of tachyonic modes leads to decoherence and spinodal 
decomposition. We compute the decoherence time, the spinodal time, the size of 
the formed domains and the homogeneous classical fields within a single domain. 



An inflationary phase in the early history of the universe can explain its present 
flatness and homogeneity as well as the anisotropy of the cosmic microwave background 

Inflation has to terminate in some non-equilibrium process which eventually leads 
to a thermal state. This is elegantly realized in hybrid models of inflation where the 
time evolution of the inflaton fleld (p controls the transition from a false vacuum state 
to the true vacuum 0. Hybrid inflation arises naturally in supersymmetric theories 
with spontaneous symmetry breaking [^. The dynamics of the inflaton fleld during the 
inflationary phase is then determined either by supergravity corrections , by quantum 
corrections |^ or by the effects of the supersymmetry breaking sector of the theory |Q. 

The phase transition at the end of inflation is of crucial importance since it determines 
the initial conditions for the further evolution towards the thermal state. Important 
quantities, which have to be obtained, are the reheating temperature and the initial 
abundances of weakly coupled particles, such as gravitinos, which are not in thermal 
equilibrium. In hybrid inflation the decay of the false vacuum takes place via spinodal 
decomposition p, similar to a second-order phase transition. Such a decay has already 
been discussed in the literature in the context of new inflation J^, and also for scalar 
0^-theory in Minkowsky space 0. 

In the following we shall study the onset of this transition in a class of hybrid 
inflation models where the evolution is controlled by the slow-roll motion of the 
inflaton fleld. This will allow us to follow in detail the breakup of the homogeneous 
inflationary phase into domains and to compute the decoherence time, the spinodal 
time, the average size of the formed domains and the homogeneous classical flelds inside 
a single domain. 

False vacuum inflation 

Hybrid inflation can be based on the superpotential |^ 

W = Wg + Ws = XT (Ml - S^) + M|(/? + ^) . (1) 



Here Wq describes the supersymmetric breaking of a global symmetry, whereas the 
Polonyi superpotential Ws accounts for supersymmetry breaking |TD|. The discrete Z2 



symmetry of the model leads to the well known domain wall problem [|Tl|. We consider 
this case only for simplicity; the extension to a continuous U{1) symmetry, where no such 
problems exist, is straightforward. The mass Ms is the scale of supersymmetry breaking, 
associated with the gravitino mass m3/2 = M|/Mpexp(2 — y^), and the constant (3 ^ 
(2 - V3)Mp, where Mp = {St^GnY^^'^ = 2.4 x 10^^ GeV is the Planck mass. Successful 
inflation requires M| -C XMq. A rotation of the flelds T and S by the small angle 



2 



^ ~ Mg/{XMQ) changes the superpotential (P to the standard form of an O'Raifeartaigh 
model 0. 

In the case of global supersymmetry one obtains from the superpotential (|1|) the 
scalar potential 



S ' 



with the ground state 



(T) = 0, |(S) 



M, 



G 



(2) 



(3) 



For |Tp > Mg/2, Vq has a local minimum at S = where it is flat in S and T. 
Supergravity corrections are included in the potential 



V 



dW ^ z*W 



dzi 



Ml 



\W\^ 



(4) 



where Zi = T, S, S", i = 1 ... 3. We choose K to be the canonical Kahler potential, 



K = which is singled out in = 2 supersymmetry [13|. Up to corrections 

0{1/Mp), the expectation values of the three fields T, S and S in the true vacuum are 



(r) = 0, |(S)| = Mg, (5) = (2 - v^)Mp . 



(5) 



The supergravity corrections distinguish between the real and the imaginary part of 
T ~ ((f + ix) I V^- For ip > (pc = Mg, V has again a local minimum at S = (cr + ip) / \/2 = 
0. The energy density is then V ~ X^Mq, which leads to an inflationary phase with 
Hubble parameter H = XMq/ (y/SMp). Here the potential is almost fiat in Lp which plays 
the role of the inflaton. The mass of p is large, rrip ~ 2X(f > 2XMg- The minimum 
with respect to S and x li^s at 5* = x = 0, but the corresponding masses are small, 
ms = 0{H) and = 0{H<p/Mp). 

In this paper we are interested in the decay of the false vacuum driving inflation, 
which takes place near the critical point (fc = Mg- This phase transition is determined by 
the dynamics of the lightest scalar fields, the inflaton (f and a, which becomes massless 
at ipc- The scalar potential for ip and a reads, up to corrections 0(1/Mp), 



V{ip, a) = aM M? 



a 



I \ 2 2 2 I 

+ X ip a + Kip 



G 



(6) 



The linear term in the potential is (9(1/Mp), with 



2^/^(2-^3) A f ^ 
M-p 
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0.6 Amg/aM^ . 



(7) 



For a certain range of parameters the time evolution of the inflaton field during 
inflation is determined by the linear term in the potential Representative values are 
Ms ~ 1.6 X 10^° GeV, which implies mg/s ^ 130 GeV, A ^ 1.3 x 10"^ and Mq ~ 10^^ GeV. 
This yields ~ 3.1 x 10^ GeV and k ~ (4.6 x lO^GeV)^. From the full supergravity 
potential (^) one obtains for the slow-roll conditions of the inflaton ip, as long as is 
sufficiently small, 

These conditions are clearly satisfied for values of -C Mp and in particular below the 
critical value ipc- The solution of the equation of motion for the homogeneous classical 
field reads, 

^,i{t) = MG-^t, (10) 

where we have fixed the origin of time by y9(0) = = Mq- The transition from the 
inflationary phase to the thermal phase starts at t = 0. 



Fluctuations near the critical point 



At <f = (fc the mass term of the cr-field changes sign. This is analogous to a second- 
order phase transition with (f{t) playing the role of a time dependent temperature of the 
thermal bath. In flat de Sitter space with metric 

ds^ = dt" - e^'^'dx^ , (11) 

the quantum fields ip{t,x) and a{t,x) can be expanded into plane waves which are mo- 
mentum eigenstates, 

^{t,x) = Vci{t) + e-i''' j j^^{a^{k)v,{t^^^^^ (12) 

a{t,x) = e-t^*/^^(a.(^)a,(t)e^*^^>at(fc)a^(t)e-^^^^) . (13) 

Note the absence of cTdit), the homogeneous part of a(t,x). This is due to the initial 
condition aci{0) = at the end of inflation. The annihilation and creation operators 
satisfy canonical commutation relations, 

K(fc),a^(fc')] = Kik),at{k')] = 5\k-k') , (14) 
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and the canonically conjugate operators to the fields and a are 

7r^(t, x) = e^^^ipit, x) , 7r„(t, x) = e^^^&{t, x) . (15) 

The assumption of a constant Hubble parameter is well justified during infiation and the 
onset of the subsequent phase transition. 

Near the critical point if = ip^ the potential for the quantum fields ip and a takes the 
form, 

Viip,a) = -^Dha' + ... . (16) 

where = AX'^k,Mg/{3H). For the parameters given above one has D = 3 AH. Note, 
that this equation also holds in other infiationary models where the infiaton is slowly 
rolling near ip = ipc, since D^t is the first term in an expansion of the cr-mass around 
t = 0. 

The mode functions ipk{t) and (Tk{t) satisfy the equations of motion, 

V^. + (k'e-'"' - ^H'^ ipk = 0, (17) 
dk + (fc^e-^^* - ^H' - DH^ a, = . (18) 

The solution of eq. (|1^ is well known from studies of density fiuctuations during inflation. 
It is given by a Hankel function 

^.W = ^(|:)'^' H«(^), (19) 

where z = kexp {—Ht)/H and 



The solution of eq. (18) has to be determined numerically. However, asymptotically, 
for t <^tk and tk, where 

k2^-2m, _ 9^2 ^ |^3^^| ^ (21) 

useful approximate analytical solutions can be foundQ. At large negative times we again 
have the Hankel function analogous to (PD|) as solution. This is the de Sitter-space initial 
condition. 



^For momenta 3H/2 < k < kmax — {D^ /{2He)y/^ ~ 3H eq. ( ^ ) has more than one solution for tk- 
We choose the smallest value for tk to construct an approximate solution. For a more detailed discussion, 
see ill. 



For t ^ tk eq. ([ISD reduces to 

dk - Dhak = 



(22) 



which describes a scalar field in Minkowski space with mass squared /i^(t) = —D^t, 



decreasing with time. The solutions are given by the Airy functions Ai and Bi [|15[. An 
approximate analytic solution of eq. (plSf ) can be obtained by matching Airy and Hankel 
functions at t = tk, 

ak{t) = CA{k) Ai{z) + CB{k) Bi{z) , (23) 
where z = Dt. For soft momenta, k < H, one finds 



For z = Dt ^ 1 the behaviour of (pk(t) is dominated by the asymptotic growth of Bi(z), 

Bi(^) ^ ^z-'/^exp(-z'/A . (25) 



% V3 

The physical meaning of such tachyonic instabilities was studied in detail by Guth and 
Pi for fixed negative mass squared, /i^(t) = —const, in the context of new infiation 0. 
More recently, the dynamics of spinodal decomposition has been investigated by several 
groups g |[. 

The situation considered here is similar to the previously discussed cases with respect 
to the occurrence of spinodal decomposition. However, a crucial difference is that the 
negative curvature of the potential is neither constant nor suddenly 'turned on'. The 
change from positive to negative curvature is, on the contrary, controlled by the slow-roll 
of the inflaton field across the critical point if = (fc- 

An important quantity is the growth of the total fiuctuation, the variance 

{a^{t)\. = {a'{t)) - {a'm , (26) 

where 

(^'W) = / 4^^e-'">k{t)\' ^ f dkhit) . (27) 



(2vr) 



Results for the spectral function fk{t) are shown in Fig. 1. The expected rapid growth 
with time is clearly visible. For small k one has fk oc /c^, since |Cb(A;)| ~ const. For 
larger momenta, k > k^ 3H, also the time t^ is large, which leads to a delayed growth 
of these modes and a rapid fall-off of fk{t) at large k. For very large momenta, k ^ k^,, 
one sees the behaviour of a massless field, i.e. fk oc k, which is redshifted with increasing 
time due to the Hubble expansion. 
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Figure 1: Spectral function fk{t) of the variance (o"^(t))y^. Thick lines correspond to 
numerical solutions; the thin line denotes an approximate analytical solution. 

The growth of the variance {a'^{t))r is shown in Fig. 2. It clearly exhibits the char- 
acteristic feature of our model, a faster than exponential growth with time. Eventually, 
this growth is terminated by the back reaction of the produced fields. The spinodal time 
at which the formation of domains is completed can be estimated by comparing the 
variance with the global minimum of the potential. 



(o-^(isp))r 



2M^ 



(28) 



From Fig. 2 one reads off t^p — 3/iJ. 

Tachyonic modes ak become classical at the decoherence time tdec- At later times the 
system is characterized by a classical probability distribution f{ak,t) for the amplitude 
(Tfc. The homogeneous mode of the system (Tci(t) = {a(t,x)) remains zero after spinodal 



decomposition ||T6 



The decoherence time t^ec is determined by the requirement that the tachyonic modes 
become classical, i.e., the product \o'k{t)7r„k{t)\ has to be larger than 1/2, the minimal 
value obtained for an oscillating mode 0. For large times, Dt ^ 1, the product is 
determined by the exponentially growing Airy function Bi(z), 



\ak{t)7T^k{t)\ ^ -D\CB{Wexp [■^{DtY^ ) ■ (29) 
For soft modes k < H, the condition \crk{tdec)'^ak{tdec) \ = Rdec ^ 1 yields the decoherence 




Figure 2: Growth of the variance (cr^(t))r with time. Full and dashed lines correspond 
to numerical and analytical results, respectively. The horizontal line denotes the true 
vacuum (o"^)r = 2M^. 



time 

tdec - 7^ (4 In {2Rdec)j . (30) 

For Rdec = 100 this implies t^ec ~3/-D~l/-f^. Results for the decoherence time obtained 
from the numerical integration are shown in Fig. 3. 



From eq. (10) and from Fig. 3 one obtains the values of the infiaton field (p and the 
variance (cr^)^/^ at the decoherence time tdec ~ 3/D, 

Mg - ^ciitdec) ~ 10^' GeV , {a\tdec))y' ~ 10^ GeV . (31) 

These values can serve as initial conditions for the further time evolution of the system 



T7| . Note, that they are different for the two fields and much larger than the frequently 



assumed initial condition H/ {2tt) 



Formation of domains 



The rapid growth of tachyonic modes leads to a breakup of the homogeneous infla- 
tionary phase into domains. This spinodal decomposition is conveniently described by 
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Figure 3: Decoherence time for tachyonic modes as function of momentum. 

means of fields averaged over a volume Vi = characterized by some smearing function 
0, [TB|. For simplicity we shall use just a momentum cutoff {ki = 2it/1), 



ai{t,x) 



e 2 



Ht 



I -— ^ {a^{k)a,{t)e^''^ + ai{k)al{t)e 

k<ki [27iy'^ ^ 



-ikx 



One then obtains for the correlation function of the smeared operator {x = \x\), 

1 



{ai{t,x)ai{t,0)) 



-3Ht 



dkk\ak{t)\ sin kx . 



(32) 



(33) 



2tT^X Jk<ki 

At time t the physical distance and the physical momentum are given hj y = x exp (Ht) 
and p = fcexp {—Ht), respectively. 

An analytic estimate for the correlation function can be obtained based on the match- 
ing of Hankel and Airy functions described above. Consider first the case of short times, 
Dt < 1, and short distances, 1/ki < x <^ 1/H. Here one finds, 

1 



{ai{t,x)ai{t,0)) 



(cos (kix) 



(34) 



This is the result expected for a massless scalar field in Minkowsky space. The wavelength 
of the oscillation is given by the ultraviolet cutoff, / = 27i /ki. 

For large times, Dt ^ 1, the correlation function is dominated by the tachyonic 
modes with k < k^. For distances x < the deviation from homogeneity is small, 

{ai{t, x)ai{t, 0)) ~ {a\t))r (l + 0{ik,xf)) . (35) 
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Figure 4: Correlation function of the smeared a-field. 

For large distances, x > l//c*, the correlation function rapidly falls off. 

Numerical results for the correlation function are shown in Fig. 4 for two different 
values of the ultraviolet cut-off ki. For a large cut-off, ki = lOH, one reads off the 
correlation length Icorr — 1/(3-?^) — l/Zc*, which corresponds to the average size of a 
domain. For cut-offs below 1/lcorr, e.g. ki = H, the expected decrease of the correlation 
function is observed. The energy density stored in spatial variations of the cr-field is 
roughly given by Pgrad ^ ^h/^lorri which is much smaller than the total energy density 
V ~ A^M^. 

So far we have ignored all back reaction effects. They will terminate the slow-roll 
motion of the infiaton field, damp the growth of the tachyonic modes and modify the 
spinodal time t^p. However, we do not expect a change of our qualitative picture with 
small decoherence time, tdec < tgp, and asymmetric classical homogeneous fields if and 
a larger than H/{27r). The determination of the properties of the final thermal state, 
in particular the reheating temperature and the energy fraction of topological defects, 
requires further work. 

Thanks to the slow-roll motion of the infiaton field across the critical point (^c we 
have been able to follow in detail the breakup of the homogeneous infiationary phase into 
domains. The tachyonic modes of the field a, which become classical before the spinodal 
decomposition ends, lead to an average domain size 0(1/H). Within a single domain 
they appear as a homogeneous classical field. 
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